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^ ■ Abstract 

The theoretical status of direct CP violation e'/e is briefly reviewed. Special 
attention is paid to the recent new consistent predictions for both the ratio 
e'/e and the AI = 1/2 rule within the standard model. In particular, two 
matching conditions resulting from the matching between the QCD and chiral 
perturbation theory(ChPT), and also some algebraic relations of chiral oper- 
■ ators are found to be very useful. It is of interest that the new predictions are 

no longer sensitive to the strange quark mass, and are also renormalization 
scale and scheme independent in the leading QCD and chiral loop approxima- 
tion with large iV c approach. The new prediction for the direct CP violation 
with the value e'/e = (20 ±4) x 10~ 4 [Im\ t /1.2 x 1CT 4 ] is consistent with the 
most recent experimental results reported by the NA48 and KTeV groups. 
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I. INTRODUCTION 

One of puzzles in particle physics is the origin and mechanism of CP violation 111 the 
standard model (SM), CP violation is described by the Kobayashi-Maskawa (KM) || phase 
in the Cabbibo-Kobayashi-Maskawa quark mixing matrix. Such a phase can arise from the 
explicit CP violation imposed in the Yukawa coupling constants of SM or originate from 
spontaneous CP violation H in some extended models of SM. One of the simplest models 
is the extension of SM with two Higgs doublets (S2HDM) Q motivated from spontaneous 
CP violation. In such an S2HDM, rich CP-violating sources can be induced from a single 
relative CP-violating phase of vacuum. These sources have been classified into four types 
Q: (i) induced single KM phase 5km ot Wolfenstein parameter rj ||; (ii) induced new 
type of flavor-dependent CP- violating phases 5f. via charged and neutral scalar interactions; 
(iii) induced super-weak type CP violation || via flavor- changing scalar interactions; (iv) 
induced neutral scalar-pseudoscalar mixing CP violation. When the source of CP violation 
is dominated by the induced KM phase, the phenomena of CP violation in S2HDM would 
be similar to the SM. 

To make a consistent prediction for the direct CP- violating parameter e'/e caused by 
the KM CP-violating phase, it is necessary to understand simultaneously the longstanding 
puzzle of the AI =1/2 rule in the kaon decays as they involve the long-distance evolution 
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of common hadronic matrix elements, where the low energy dynamics of QCD shall play 
a crucial role for a consistent analysis. During the past few years, both theoretical and 
experimental efforts on direct CP violation in the kaon decays have been made important 
progresses. On the experimental side, two improved new experiments [|7|,§ with higher pre- 
cision have reported results which are consistent each other at the 1 — a level. On the 
theoretical side, several groups P-fZll have made detailed calculations for the ratio e'/e. Re- 
cently, there have been some interesting developments || which are mainly based on QCD 
of quarks and chiral perturbation theory (ChPT) p2| . |23| at low energies for mesons. Conse- 
quently, it has reached an agreement between the experimental results and the theoretical 
predictions. 

To reach a consistent theoretical prediction for the direct CP-violating parameter e'/e, it 
must satisfy, at least, the following simple criteria: (i) the prediction should be renormaliza- 
tion scale independent; (ii) the AI =1/2 rule can be reproduced. To satisfy the first criteria 
(i), one should be able to solve the matching problem between perturbative QCD used to 
calculate the Wilson coefficient functions c,(/i) and effective theories applied to evaluate 
the hadronic matrix elements of operators < \Qi\ > (//). To reach the second one (ii), the 
effective theory should well describe the low energy dynamics of QCD in weak kaon decays. 
One of such attractive effective theories is the chiral perturbation theory (ChPT) inspired 
from 1/N C expansion [23,25], which has actually provided a successful description on many 
processes [p3fl . In ref. |J, we have further shown that the ChPT with a functional cutoff 
momentum scheme can lead to a consistent prediction for the direct CP- violating parameter 
e'/e and the AI = 1/2 rule. 



II. BRIEF REVIEW ON THEORETICAL STATUS 

The present status of theory versus experiment is summarized as follows: 
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It is seen that great theoretical efforts have been made by many groups as a consistent 
prediction for e'/e and A J = 1/2 rule play an important role for understanding the origin 
and mechanism of CP violation as well as the low energy dynamics of QCD. Where different 
approaches have been adopted to obtain the above theoretical results: 

Beijing's : The ChPT inspired from 1/N C expansion has been adopted to calculate the 
chiral loop contributions. The claculating scheme is based on the one with cutoff momentum 
M suggested by Bardeen, Buras and Gerard [pEfl , but it is generalized by taking the cutoff 
momentum as the function of the renormalization scale /i, i.e., M = M(/x), instead of naively 
identifying the cut-off momentum to the QCD running scale \x. It is such a simple extention 
that enables us to make the new precitions for both the direct CP-violating parameter 
e'/e and the AI = 1/2 rule be renormalization scale and scheme independent as well as be 
consistent with the experimental data The new predictions also nicely improved our previous 
ones presented in ref. |[27|| . 

Dortmund's [It]]: Chiral loops are calculated by separating the factorized and non- 



factorized diagrams |f28|| . The AI =1/2 rule may be reproduced ||29|| . While the matching 



procedure needs to be improved as the presence of the quadratic cutoff may induce a match- 
ing scale instability. 

Dubna's Chiral loops are regularized via the heat-kernel method in the ENJL 

framework up to 0(p 6 ). The renormalization scheme dependence needs to be improved. A 
phenomenological fit of the AI = 1/2 rule also results in a big uncertainty on the matrix 
element (Qe). 



Granada-Lund's ||12|| : The X-boson method within the ENJL model is used to improve 
the matching between long- and short- distance components. The AI = 1/2 rule can be 
reproduced. Both renormalization scale and scheme dependences have been improved. While 
Non-FSI chiral corrections need to be further considered. 



Montiepeller's ||13|| : Effects to the matrix element (Qq) from the qq component of the 
scalar meson was considered by using QCD sum rule. The question is how to explicitly 
separate the qq from gluon component of the scalar meson. 



Munich's ||14|| : A phenomenological 1/N C approach is adopted, where some of the 
matrix elements were obtained by fitting the AI = 1/2 rule at fi — m c = 1.3 GeV. The 
matrix elements (Qq) and (Q$)2 relevant to the direct CP- violating parameter e'/e remain 
undetermined and are taken around their leading 1/N C values. 



Roma's [|15|| : Present lattice results remain unreliable as large renormalization uncer- 
tainties at the matching scale between the lattice and continuum results. The present lattice 
calculations can only use the lowest order chiral perturbation theory to evaluate the K — > tttt 
amplitude. The recent result for e'/e was estimated by taking Bq to be the VSA result varied 
by a 100% error. 

Taipei's fll6| : Effective Hamiltonian approach was adopted in conjunction with general- 
ized factorization for hadronic matrix elements. The non-factorizable effects were considered 
by introducing phenomenological parameters. Some assumptions are needed to fix the phe- 
nomenological parameters. 

Trieste's ||17|| : The chiral quark model has been used to evaluate the hadronic matrix 
elements which were assumed to be matched to the NLO short-distance Wilson coefficients 
at \i = 0.8 GeV. The AI = 1/2 rule can be reproduced at that point. The problem of 
renormalization scale and scheme dependence needs to be improved. 
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Valencia's [|18|| : A dispersive analysis a la Omnes is used to obtain the FSI effects 
relative to the leading 1/N C amplitudes. The FSI effects alone cannot reproduce the AI = 
1/2 rule. More comments on this approach may be found in ref. [[J(J . 

In most of the approaches, the matching procedure remains the main problem. This is 
because the energy scale M of long-distance operator evolution from meson loops must in 
general be smaller than the chiral symmetry breaking scale Af, i.e., M < Af ~ 1 GeV, while 
the energy scale \x of the short-distance operator evolution from perturbative QCD should 
be above the confining scale, i.e., \i > 1 GeV. Naively identifying the scale M in ChPT 
to the scale \i in perturbative QCD may become inappropriate. In the following, we will 
focus our discussions on the recent new consistent predictions for the direct CP-violating 
parameter e'/e and the A J =1/2 rule presented in ref. [f|. 



III. LOW ENERGY DYNAMICS OF QCD 

To treat the low energy dynamics of QCD, our starting points are mainly based on the 
following basic considerations: 

• In the large N c limit but with the combination a s N c = a being held fixed. The QCD 
loop corrections which are proportional to a s are then corresponding to a large N c 
expansion, a s ~ l/N c |2^ 



Chiral flavor symmetry SU(3)l x SU(3)r is supposed to be broken dynamically due 
to attractive gauge interactions, namely the chiral condensates < qq > exist and lead 
to the Goldstone-like pseudoscalar mesons n, K, rj. The chiral symmetry breaking 



scale A j is characterized by the condensate, Af ~ 4iT^—2<qq> jr ~ 1 GeV with 
r = m^ o /m (m = (m u + m^)/2). 

The low energy dynamics of QCD in large N c limit is considered to be described by 
the chiral Lagrangian. The ChPT is going to be treated as a cut-off effective field 
theory which may be regarded as a consistent theory in a more general sense BT] . 



The cut-off momentum M is expected to be below the chiral symmetry breaking scale 
A/- 

The chiral meson loop contributions are characterized by the powers of p 2 /A 2 with 
Af = Airf . Here f 2 ~ —2 < qq > jr ~ N c is at the leading N c order and fixed by the 
7r decay coupling constant / ~ F^. Thus the chiral meson loop contributions are also 
corresponding to a large N c expansion of QCD, p 2 /A"j ~ 1/N C ~ a s . Therefore both 
chiral loop and QCD loop contributions must be matched each other, at least in the 
sense of large N c limit. Thus the final physical results should be independent of the 
renormalization scale and the calculating schemes. 

The cut-off momentum M of loop integrals is in general taken to be a function of /i, 
i.e., M = M(/i) 0, instead of naively identifying it to the renormalization scale \x 
appearing in the perturbative QCD in large N c limit. M = M(/i) may be regarded as 
a functional cut-off momentum, its form is determined by the matching between the 
Wilson coefficients of QCD and hadronic matrix elements evaluated via ChPT. 
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In a word, we are going to treat the ChPT with functional cut-off momentum M(/i) 
scheme, as a low energy effective field theory of QCD in the large N c approach. 



IV. CHIRAL REPRESENTATIONS AND RELATIONS OF OPERATORS 

The general form of the chiral Lagrangian is constructed based on the chiral flavor 
symmetry SU(3)l x SU(3)r and can be expressed in terms of the expansion of the momentum 
p and the quark mass m q to the energy scale A x = 0(1) GeV at which nonperturbative 
effects start to play the rule. For our purpose, here we only use the chiral Lagrangian which 



is relevant to the K — ► rnr decays (for the most general one, see ref. [23] ) 



C eff = J —{ triD^D^U) + j^tr{ln U ] - In U) 2 + r tr{MU ] + UM ] ) 



+r||tr (p^D»U{M^U + Wm) (1) 



X5_ 

A 2 

x 

+r 2 ^-tr (M ] UM ] U + MU ] MU ] ) + r 2 ^-tr{M^M) } 
A x ' A x 

with DpP = dfJJ — ir^U + iUl^, and M. = diag(m u , nid, ra s ). and r M are left- and right- 
handed gauge fields, respectively. The unitary matrix U is a non-linear representation of 
the pseudoscalar meson nonet given as U — e tYl ^ with IT = 7r a A a and tr(X a \b) = 25 a b- 
Here we keep the leading terms at large N c limit except the anomaly term which arises 
from the order of 1/N C . Note that in order to make clear for two independent expansions, 
namely 1/N C expansion characterized by p 2 /A. 2 in the large N c limit, and the momentum 
expansion described by p 2 /A 2 i , it is useful to introduce the scaling factor A x ~ 1 GeV and 



to redefine the low energy coupling constants Lj introduced in ref. |23| via Lj = x% f 2 
and Hj = Kj / 2 /4A 2 , so that the coupling constants Xi (i — 3, 5, 8) and A x are constants in 
the large N c limit and the whole Lagrangian is multiplied by f 2 and is of order N c except 
the U(l) anomalous term. Numerically, one sees that Xi = 0(1) for A x = 1 GeV. 

with the above chiral Lagrangian, the quark currents and densities can be represented 
by the chiral fields 

ST f 2 



X5_ 

2A 2 



dHf*M - M ] d^U + d^UM^U - U^MU^d^u) } ij , (2) 



(3) 



Similarly one can obtain the right-handed currents and densities. With these definitions, 
all the current x current and density x density four quark operators can be reexpressed in 
terms of the chiral fields, we may call such chiral representations of four quark operators Qi 
as chiral operators denoted by Qf correspondingly. 

In the standard model, it has been found that the AS = 1 low energy (/x < m c ) effective 
Hamiltonian for calculating K —>■ titt decay amplitudes can be written as 
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n^ 1 = ^Kfjk{ J i)Q i { l i) , (/i<m c ) (4) 

with Qi the four quark operators. Note that only seven operators are independent as the 
linear relation Q 4 = Q 2 — Qi + Q3. Q7 and Q$ originate from electroweak penguin diagrams. 
Ci(/i) are Wilson coefficient functions Cj(/i) = Zi(ji) + ryi(fi). Where r = —X t /X u with 
X q = V* s V q d- The hard task is for calculating the hadronic matrix elements (Qi(f))i for 
< A x = 1 GeV which is at the order of chiral symmetry breaking scale. This is because 
perturbative QCD becomes unreliable in such a low energy scale. It has been shown || 
that the ChPT with functional cut-off momentum could provide a powerful way to evaluate 
(Qi(n))i when fi < A x . The procedure is as follows: Firstly one represents the current 
x current or density x density four quark operators Qi by bosonized chiral fields from the 
chiral Lagarangian, i.e., Qf, then calculate loop contributions by using the functional cut-off 
momentum scheme. Finally, one matches the two results obtained from QCD and ChPT 
with functional cut-off momentum by requiring scale independence of the physical results. 
Qf can be written as the following form |9j 

qx + Fc . = _f fr (XzrfdpU) tr (A (1) f/W) + 0(1/AJ) , 
Qf + H.c. = -f tr (x & U^d^UX {l) U^d»u) + 0(1/Aj) , 
qx + E c _ = _f tr (x^d^U) tr (U^U) + 0(1/ A 2 X ) , 
qx + Ha = _ f 4 tr (x^rfd^u) + 0(1/AJ) , 

qx + Ha = _f tr (Xetfd^u) tr (Ud*rf) + 0(1/ A 2 X ) , (5) 

qx + F c . = +/ 4 tr ^d^u) + 0(1/ A x ) , 

qx + H _ c _ = J_qx _ 3 /4 fr (Xstfd^u) tr (xPUPrf) + 0(1/ A 2 X ) , 

Q x s + H.c. = -\q\ + fAr (yji\0)u) 

+ / 4 r 2 |||trA 6 (rf\V>Ud li rfd tl U + d^U ] d»UU ] X m u) 

x 

+ A 2 -||2rtrA 6 (U ] X {1) U M^U + U ] MU ] X {1) U) +0(1/Aj). 

x 

with the matrix A^^diag. (1,0,0). Thus loop contributions of the chiral operators Qf can 
be systematically calculated by using ChPT with functional cut-off momentum. 

For K — > 7T7T decay amplitudes and direct CP- violating parameter e'/e, the most im- 
portant chiral operators are Qf , Qf, Qf and Qf. In fact, the chiral operators Qf and Qf 
decouples from the loop evaluations at the p 2 order, i.e. 

Qf = Qf = Q (6) 

which can explicitly be seen from the above chiral representations due to the traceless factor 
tr = when ignoring the singlet U(l) nonet term which is irrelevant to the Kaon 
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decays. Here Ud^W = A a ^X a may be regarded as a pure gauge. As a consequence, it implies 
that at the lowest order of p 2 , we arrive at two additional algebraic chiral relations 

qx = qx_ q x = _ /4 tr ^d^&^U) + 0(1/ A 2 X ) . (7) 

and 

Q X e = ~ (^) (Q X 2 ~ Qi) = (^) f tr fad^tru) (8) 

Notice that the mass parameter r is at the same order of the energy scale A x , and x$ is a t 
order of unit, thus the leading non-zero contribution of Qq is at the same order of Q\ and 

Ql 



The above algebraic chiral relations mentioned also as Wu-relations in the literature |32 



were first derived in ref. [p7[] , they have been checked from an explicit calculation up to 



the chiral one-loop level by using the usual cut-off regularization ||26|| . It was based on this 



observation that the ratio e'/e was predicted to be large enough ( e'/e — (1—3) x 10 

) for observation. In ref. it has further been shown that the algebraic Wu-relations of 
the chiral operators survive from loop corrections and can provide a consistent prediction 
for both the direct CP- violating parameter e'/e and the AI = 1/2 rule when applying for 
an appropriate matching procedure between QCD and ChPT. 



V. MATCHING APPROACH 

The short-distance evolution of quark operators is performed from perturbative QCD. 
When the energy scale fi is high, m w > jj, > m b , there are eleven independent operators Qi 
{i = 1, • • • , 11). When the energy scale \i runs down to below the bottom quark mass m b 
and above the charm quark mass m c , i.e., nib > ^ > m c , the operator Qu decouples and 
operator Q w is given by the linear combination Q w = —2Q\ + 2Q 2 + Qz — Q\. Once the 
energy scale fi goes down to below m c but above the confining scale or the energy scale A x , 
i.e., m c > /i > A x , two operators Q$ and Q4 become no longer independent and are given 
by the linear combination Q 9 = Q 2 + Q\ and Q 4 = Q 3 + Q 2 — Q\. Thus there are only seven 
independent operators below m c and above A x . The one- loop QCD corrections of the quark 
operators at the energy scale just above the energy scale A x is given by 

2 

Q^q) = QM - 3^ ln(£f) QM + 0(1/JV C ) , (9) 

2 

Q 2 (^)=Q 2 (/i)-3^1n(^f)Q 1 ( /U ) 
Air fi^ 

- \ ^ in(4) QM - 1 r ln (4) QM + 0{1/N e ) , (10) 

3 47T \l z 3 47T \i l 

2 



Art /j 

^ M 4) QlW _^i n( ^ 

47T V /i V ' 47T v fJL 



^ ln(^f ) Q 4 (/i) - t4 I*® QM + 0(1/N C ) , (11) 
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2 

Q*fa) = QM-T L H ! %)QM 



,2 



+ [3(N C - l/N c ) - 1]^ ln(^) g 6 (/i) + 0(1/JV C ) , (12) 



,2 



g 8 (/i Q ) = QM + [3(iV c - l/JV c ) - 1]^ ln(^) Qs(fi) , (13) 



and 



Qsfa) = QM ~ T T- ln( 4) <W -\%- In(4) QeM + 0(l/iV c ) , (14) 

2 

Q 5 (/i Q ) = g 5 (/i) + 3^ ln(^f ) g 6 (/i) + 0(l/iV c ) , (15) 

2 

g 7 (/i Q ) = g 7 (/x) + 3^ in(^f ) g 8 (/i) + o(i/tv c ) . (ie) 

47T /i z 

It is then clear that 

i) In the large N c limit, Q±, Q2, Qa, and Qq form a complete set of operators under QCD 
corrections. 

ii) The evolution of Q 8 is independent of other operators and only caused by loop cor- 
rections of the density. 

iii) The operator Q 3 is given by the linear combination Q 3 = Q A — (Q 2 — Qi). The 
operator Q 5 is driven by the operator Qq, and the operator Q 7 is driven by the operator Q 8 . 

When the energy scale fi approaches to the confining scale, or /i < A x ~ Ap ~ lGeV, 
as we have discussed in the above sections, long-distance effects have to be considered. The 
evolution of the operators Qi(fi) when jj, < A x is supposed to be carried out by the one of 
the chiral operators Q}{M(n)) in the framework of the functional cut-off ChPT truncated 
to the pseudoscalars. To be treated at the same approximations made in the short- distance 
operator evolution of QCD, we should only keep the leading terms (i.e., quadratic terms of 
functional cut-off momentum) and take the chiral limit, i.e., m^rn^ << A 2 F . The evolution 
of the operators Q\ and Q\ is simply given by 

- Q?(MM) = Qf(MM) - 2 ( M2 M ~ M V)) qx (m(//)) > (17) 



A 



+ M M -MW m _ Q})(M{p , )h (lg) 

where A F = AttF = 1.16 GeV with F the renormalized one of /. The operators Qf (i = 
4, 6, 8) can be written as follows 

Qa(jm) - Qi(M^)) = (Q* - g?)(M(/i)) , (19) 

QM -> g*(/i, M(jm)) =(l + 3(N c - 1/AQ ^ In(4)l OSfe M M) , (20) 



47T /i 



g 8 (/i) - C%(jm, M(jm)) =(l + 3(N c - 1/AQ g ln(^)) Q*^, M(jjl)) . (21) 
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where the explicit /i-dependence of the operators Qq(/i, M(fi)) and Qs(n, M(fj)) arise from 
the running quark mass and behavors like l/(m s (/i) + m(/i)) 2 . Notice that the operators 
and Q§ decouple from the evolution, namely Q% = and Q* = 0. Thus the independent 
operators are reduced once more in the long-distance operator evolution when fi < A x due 
to the algebraic chiral operator relations. 

Matching the loop results evaluated from QCD with the ones from the ChPT with 
functional cut-off momentum at the energy scale A x with large N c expansion and requiring 
/i-independence in the large N c limit, i.e., ^-Qi(//q) = 0, we arrive at the two matching 
conditions ||, 

= 9^' 22 

C%fa,M(ri) = ~((%-Q?)(M(ri), /i<A x (23) 
When combining the second matching condition with the chiral Wu-relation 

QU» x ,M(ri)^(-^yQX-QX)(M(fi)), /i<A x (24) 
R x = R(jx ~ fi x ) ~ ml/m(n x ) - ^ 2 K /(m s + m)(/x x ) , (25) 
we are able to fix the strange quark mass M 




Rjxl _ ii 

A x - 2 



m s (ji x ) ~ 196MeV , (26) 



Note that the coupling constants X5 an d r must be replaced by the corresponding renormal- 
ized ones xl an d -R(a*) as their loop corrections are at the subleading order. Here we have 
also used the result A x = 1.03 V / X5 GeV which is fixed from the ratio of the kaon and pion 
decay constants. 

After integrating the first matching condition, we find that the /i-dependence of the 
functional cut-off momentum M(//) can be written as || 

A 2 F Air 47r /io 

which shows that after imposing the matching condition for the anomalous dimensions be- 
tween quark operators Qi(fi) in QCD and the corresponding chiral operators Qf(M(fj,)) 
in the ChPT with functional cut-off momentum, the dimensionless ratio M 2 /A 2 F is only 
related to the strong coupling constant a s and becomes scheme-independent, which im- 
plies that the long-distance operator evolution in ChPT with functional cut-off momentum 
can be carried out by using any approach. Where fio is the low energy scale arised as 
the integrating constant. In general, we have fi > Aq C d- To fix the value of /i , we use 



M ~ yU . Thus /io (or a s (/xo) ) is determined via /^o — AF^3a s (fj, )/8ir. Using the definition 
a s (/x) = 67r/[(33 — 2nf) ln(/i 2 /Ag CD )] with rij = 3, fi is found to be (for Aq C d = 325 ± 80 
MeV) 
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Ho ~ 435 ± 70MeV or a s (/j )/27r ~ 0.19±g;g| . (28) 

Thus the functional cut-off momentum M(/i) at /i = A x yields the following corresponding 
value 

M x = M(/i = A x ~ lGeV) ~ OJll^GeV . (29) 

which provides the possible allowed range of the energy scale where the ChPT with functional 
cut-off momentum can be used to describe the low energy behavior of QCD at large N c limit. 



VI. LONG-DISTANCE EVOLUTION OF CHIRAL OPERATORS 

We may now adopt the matching conditions and algebraic Wu-relations of the chiral 
operators to investigate the long-distance operator evolution. The AS = 1 low energy 
(/i < A x ) effective Hamiltonian for calculating K — > ixix decay amplitudes may be expressed 
as follows 

Htf f =1 = ^K £ c t (A x )Qf(M(A x )) (30) 



=1,2,4,6,8 



It is seen that the first matching condition enables us to sum over all the leading terms 
via renormalization group equation down to the energy scale /i , and the second matching 
condition together with the algebraic chiral operator relations allows us to evaluate the 
penguin operators Q%(M) and Qq(M) from the operators Q\(M) and Q^M). So that 
the operators Q*(M) and Q*(M) form a complete set for the operator evolution below the 
energy scale /i ~ A x ~ 1 GeV, or correspondingly, below the functional cut-off momentum 
M(/jl ~ A x ) ~ 0.7lto;j2 GeV for Aq C d = 325 ± 80 MeV. By choosing a new operator basis 
Q±(M((i)) = Q2(M(fi)) ± Q2(M(fJ.)) with the anomalous dimension matrix for the basis 

(Q x -,Ql) 

a s ( -9/2 0\ 
7= 27(-3/2 3j ^ 

and following the standard procedure of the renormalization group evolution with the initial 
conditions for the Wilson coefficient functions: c_(A x ) = c 2 (A x ) — c x (A x ) and c + (A x ) = 
C2(A X ) + ci(A x ) , we find in the leading logarithmic approximation that 

Q*(M(A X )) =r 1 - l/2 ri-{M x )Q^) , (32) 

Qim^))=v^m(M x )Ql(Ho) + l(^ 1/2 -v x /3 ) m(M x )Ql(Ho) , (33) 

with i] x = a s (A x )/a s (ii ), and 

Q^M = g^(o) + ^^g^(o), (34) 

QIM = Ql(0) - ^ g* +( 0) + ^ Q-(0) , (35) 
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where rj-(M x ), r]i(M x ) and r]2{M x ) represent the finite meson mass contributions 



3 < - X , 

Al [ + m 2 



7y_(M x )~l+ 4 + 



g, 1 + i^!n( 1 + ^ 



, 2W , 1 + <#4lnfl + ^)V (36) 



Ml 



Numerically, we use m ~ 300MeV, m K = 0.495 GeV and = 0.137GeV. When the QCD 
scale takes the value Aqcd = 325 ± 80 MeV with the corresponding low energy cut-off 
momentum fj,Q ~ 435 ± 70 MeV, we have || 

Q*(M(A X )) = (3.1718S) Q*(0) = Q 4 X (M(A X )) , (37) 

Q^(M(A X )) = (0.55S") <&(0) + (0.81SS) <£(0) , (38) 

Q*(/i x , M(A X )) = -Hgx (M(Ax)) = _(i 7 .44l|f 7 ) Q* (0) , (39) 
33 A 2 

M(A X )) = - * (Ql + Q*)(0) = 19.18 (<# + Q*)(0) . (40) 

with these analyzes, we are able to present our numerical predictions for the isospin ampli- 
tudes and the direct CP- violating parameter e'/e. 

VII. NEW PREDICTIONS FOR e'/e AND AI = 1/2 RULE 

The direct CP violation e'/e in kaon decays arises from the nonzero relative phase of 
isospin amplitudes A and A 2 



uj I ' ImAi ImA 



(41) 



e y/2\e\ \ReA 2 ReA . 
with uj = ReA 2 / ReA^ = 1/22.2, we arrive at the following general expression 

7 = if W^ /mA << ft °-^ (42) 

where A/ are the K — > tttt decay amplitudes with isospin / with 

A I cos6 I = (nn\H? f s f\K) = ^K £ Q(A x )i?e(Q*(M(A x ))> 7 (43) 

i=l,2,4,6,8 

The CP-conserving amplitudes are given by 

ReA cos6 = % Re\ u « i (A x )ife(Q?(M(A x ))) 

v2 i=l,2,4,6,8 

~ ^ iteA^*_(A x )ife(g* (M(A X ))> + iz + (A x )i?e(Q* (M(A X ))> 
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+ z 4 (A x )i?e(Q*(M(A x ))> + ^ 6 (A x )i?e(Q*(M(A x ))> ] , (44) 
ReA 2 cos5 2 = ^ ReX u £ z, l {A x )Re(Q}{M{A x ))) 2 (45) 

v2 i=l,2,8 

" % Re\ u [^z4A x )Re(Ql(M(A x ))) 2 + \z + {A x )Re(Ql{M{A x ))) 2 ] , 

and the hadronic matrix elements concerned for calculating the ratio e' / e are 
ho = (cos^o)- 1 £ yi (A x )Re(Qf (M(A X ))) (1 - £l IB ) 

1=1,2,4,6,8 

~ (cos5 )- 1 2/ 6 (A x )i?e(Q^ (M(A X ))) (1 - Q IB ) , (46) 

i=l,2,8 

~ (cos<5 2 )- 1 2/8 (A x ) J Re(g^ (M(A X ))) 2 , (47) 
where we have taken into account the possible isospin breaking effect Qj B with the most 



recent refined calculation in [34]: Qjb — 0.16 ± 0.03, which is smaller than the previously 



estimated value VLjb — 0.25 ± 0.1 [[35]], but with a large error [fffi] . The CKM factors -ReA 



and ImXt are given in the Wolfenstein parameterization E] as follows 

ReX u = Re(V: s V ud ) = X , ImX t = Im(V*V ta ) = A 2 X 5 V (48) 

To obtain the numerical results, we use the following reliable values for all relevant 
parameters: A QCD = 325 ± 80 MeV, /x = 435 ± 70MeV, A x = l.OGeV and A F = 1.16 
GeV; z_(A x ) = (z 2 - Zl )(A x ) = 2.181^;^, z + (A x ) = (z 2 + Zl )(A x ) = 0.685 =F 0.029, 
Z4 (A X ) = -(0.012±0.003) and z 6 (A x ) = -(0.013±0.003), as well as y 6 (A x ) = - (0.113lffi) 

and y 8 (A x )/a = 0.158lg;g|§; (Q-(0))o = 36.9 x 10 6 MeV 3 , (Ql(0)) o = 12.3 x 10 6 MeV 3 , 
(Ql(0)) 2 = 34.8 x 10 6 MeV 3 and (Q%(A X , 0)) 2 = 328.8 x 10 6 MeV 3 . Note that for the Wilson 
coefficient functions, we only use the leading order results at one-loop level for a consistent 
analysis since the chiral operators have only been evaluated up to the leading order at the 
chiral one- loop level, namely at the order of 1/N C ~ M 2 / A 2 F ~ a s in the large A^ c approach. 
Their values can be read following the calculations in refs. [0. The numerical values at 
fj, = A x are regarded as the 'initial conditions' for the chiral operator evolution and read 
for Aqcd — 325 ± 80 MeV. For the hadronic matrix elements of chiral operators at cut-off 
momentum M = take their values at the tree-level. For the CKM matrix elements, there 
remain big uncertainties arising from the single CP-violating phase, two matrix elements V u b 
and V c b, or the corresponding Wolfenstein parameters r], p and A. For a numerical estimate, 
we take ReX u = 0.22 and ImX t = 1.2 x 10~ 4 as the central values | 38| . With these input 
values, we obtain the isospin amplitudes 

ReA = (2.56lg;P) x 10" 4 (cos c^ 1 MeV = (3.10+^ 4 ) x 10" 4 MeV , (49) 
ReA 2 = (0.12 =p 0.02) x 10" 4 (cos^)" 1 MeV = (0.12 t 0.02) x 10" 4 MeV , (50) 

which is consistent with the experimental data: ReA = 3.33 x 10~ 4 MeV and ReA 2 = 0.15 x 
10" 4 MeV . Here the final state interaction phases, 5 = (34.2 ±2.2)° and 6 2 = (-6.9 ±0.2)° 
[ 39| have been used. Simultaneously, it leads to a consistent prediction for the direct CP- 



violating parameter e'/e 
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From the above analyzes, it is noticed that 

1. The main uncertainties for the predictions arise from the QCD scale Aqcd (or the 
low energy scale /io) and the combined CKM factor Im\ t . Nevertheless, the uncertainties 
from the energy scale may be reduced from comparing the predicted isospin amplitudes Aq 
and A 2 with the well measured ones. It is seen that the results corresponding to the large 
values of Aqcd > 325 MeV appear not favorable. 

2. Considering from the isospin amplitude A 2 , the ratio e'/e favors the low values 



~ 16 x 1(T 4 (ImX t /1.2 x 1(T 4 ) (52) 



e 
e 

while from the isospin amplitude Aq, it favors the high values 



- ~ 24 x 1(T 4 (ImX t /1.2 x 10" 4 ) 



(53) 



From the ratio of two amplitudes ReA^/ ReA 2 , i.e., the AI = 1/2 rule, the ratio e'/e favors 
the middle values 

£ > 

- ~ 20 x 1(T 4 (ImX t /1.2 x 10~ 4 ) (54) 

3. The above results are renormalization scheme independent as the consistent matching 
between QCD and ChPT occurs at the leading one- loop order of l/N c ~ a s ~ 1/Af< around 
the scale A x . The renormalization scheme dependence arises from the next-to- leading order 



of perturbative QCD which could become substantial for some of the Wilson coefficient 
functions when the renormalization scale /x runs down to around the scale A x = lGeV. 
For the long-distance evolution, the scheme is fixed by the ChPT with functional cut-off 
momentum. For matching to this scheme, it is useful to introduce a scheme independent 
basis for the perturbative QCD calculation of short-distance physics. Then applying our 
above procedure to find out the matching conditions at the next-to- leading order 1/N% ~ 
a 2 s ~ l/A 4 ?. To work out the scheme independent basis in QCD, it may be helpful to adopt 



the method discussed in ref. 41 and use the cut-off momentum basis. 



In summary, our new consistent prediction for the direct CP- violating parameter e'/e is 
j = (20 ± 4) x 10- 4 ( 12 J ^q_ 4 ) = (20 ± 4 ± 5) x 10" 4 (55) 

where the first error from the low energy scale and the second one from the CKM factor 
ImXt- Our prediction is consistent with the most recent results reported by the NA48 
collaboration at CERN and the KTeV collaboration at Fermilab || 

Re(e'/e) = (20.7 ± 2.8) x 10~ 4 (2001 KTeV) [8] (56) 
Re(e'/e) = (15.3 ± 2.6) x 10~ 4 , (2001 NAA8) fTfl ; (57) 

as well as with the world average 

Re(e'/e) = (17.2 ± 1.8) x 10~ 4 (World Average 2001) (58) 
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VIII. CONCLUSIONS 



Let me briefly summarize the new predictions by Beijing group: 

• We have Clearly made a bi-expansion: 1/N C ~ a s ~ M 2 /Aj and p 2 / h? x ( mjj/ 'A|) by 
introducing the iV c -independent energy scale. The leading non-zero contributions for 
relevant chiral operators are all at the same order of \fW c - 

• Chiral representation of operators allows us to get some useful algebraic chiral relations 
which reduces the independent operators and relates the ratio e'/e to the AI =1/2 
rule. 

• Chiral loops have been calculated by using the functional cutoff momentum (FCOM) 
(M(//)) scheme. The form of M(fi) is determined by matching between short- and 
long-distance contributions. 

• The leading order oil/N c between chiral one loops and QCD one loops has been found 
to be well matched at large-A^ limit. In this sense, ChPT does well describe the low 
energy dynamics of QCD. 

• Two useful matching conditions have been resulted, which makes the predictions for 
the ratio e'/e and the AI =1/2 rule are insensitive to the strange quark mass. 

• The results at the leading order of loops are both renormalization scale and renormal- 
ization scheme independent. 

The present applications of ChPT also indicate that the ChPT with chiral Lagrangian 
obtained based on the chiral flavor symmetry SU(3)l X SU(3)r and inspired from 1/N C 
approach is a consistent and useful effective theory. With the functional cutoff momen- 
tum scheme, the ChPT may well describe the low energy dynamics of QCD. We expect 
that the theoretical uncertainties can be further improved within the framework of ChPT. 



Nonetheless to say, the flavor symmetry has plaid an important role on flavor physics [42 . 
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